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In this work we perform atomistic model potential molecular dynamics simulations 
by means of state-of-tlie art force-fields to study the implantation of a single Au 
nanocluster on a Polydimethylsiloxane substrate. All the simulations have ben per- 
formed on realistic substrate models containing up to ~ 4.6 millions of atoms having 
depths up to ~90 nm and lateral dimensions up to ~25 nm. We consider both 
entangled-melt and cross-lin ed Polydimethylsiloxane amorphous structures. 

We show that even a single cluster impact on the Polydimethylsiloxane substrate 
remarkably changes the polymer local temperature and pressure. Moreover we ob- 
serve the presence of craters created on the polymer surface having lateral dimensions 
comparable to the cluster radius and depths strongly dependent on the implantation 
energy. 

Present simulations suggest that the substrate morphology is largely affected by 
the cluster impact and that most-likely such modifications favor the the penetration 
of the next impinging clusters. 
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I. INTRODUCTION 



Electronic devices fabricated on bendable plastic substrates are widely believed to have 
great potential for applications in electronics. Much progress has been made in this area, 
in particular in the fabrication of devices and circuits, ranging from electronic paper-like 
display devices^"- to sensor skins^, circuits suitable for radio frequency identification tags^>^, 
no n- volatile memories?, photovoltaics^i^ , smart clothing-^, and actuators^^iii. 

Stretchable electronics represents a much more challenging class of systema^^ with ap- 
plications mainly in the biomedical area, where circuits must be integrated with biological 
tissues. Stretchability demands that the circuits have the capacity to absorb large strain 
without degradation in their electronic properties. Materials used in stretchable electronics 
range from entropy elastic elastomers to energy elastic solids. Stretchable circuits are made 
of diverse materials that span more than 12 orders of magnitude in Young elastic modulus. 

In this work we focus, in particular, on stretchable electrodes, formed by metallic 
paths inserted into elastomeric substrate, mainly fabricated via the direct metallization 
of elastomers^^. Polydimethylsiloxane (PDMS) is the most used elastomer substrate be- 
cause it combines biocompatibility with suitable mechanical properties'^. Moreover, PDMS 
can be easily synthesized by available laboratory techniques and it is not attacked by 
most process chemicals^. Amorphous PDMS can be synthesized in different microscopic 
structures, depending on the applications, including entangled-melt (EM) films (composed 
by a network of chains strongly wrapped around each other) as well as cross-linked (CL) 
structures in which individual chains are chemically bonded to each other through suitable 
cross-linker moleculesA^ . 

The metallization of PDMS to produce micrometric conductive pathways is typically 
obtained by metal vapor depositiort^ or by metal ion implantation' One of the main 
drawbacks of the metal vapor deposition technique is the delamination of the conducting 
layers even at very low deformations. This problem is solved by metal ion implantation, 
which provides improved adhesion between the conducting layers and the polymeric sub- 
strate. In this technique nobel metal ions are implanted with energies in the range of the 
KeV/atomii»^, thus forming a conductive layer just below the substrate surface. However, 
it has been shown that upon the ion implantation the insulating polymeric substrate is likely 
degraded by carbonization and the presence of charges^. 
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Recently a novel implantation technique named " Supersonic cluster beam implantation" 
(SCBI) has been proposed^^i^. SCBI consists in pointing a collimated beam of neutral 
metallic clusters towards a polymeric substrate. The size distribution of the clusters ranges 
from 3 nm to 10 nm and the typical implantation kinetic energy is about 1.0 eV/atom. Even 
if the kinetic energy comparatively lower than in the ion implantation technique, the neutral 
clusters are able to penetrate up to tens of nm into the polymeric target forming a buried 
conducting layer and avoiding charging and carbonization. The mechanism driving SCBI 
is the so-called " clearing-the-way" effect^""—, namely the front atoms in the cluster collide 
with target atoms conveying them sufficient momentum to clear the way for the following 
cluster atoms. As a consequence, the low-energy cluster penetration depth is increased with 
respect to the high-energy single ion penetration. The "clearing the way" effect has been 
confirmed and formalized by classical molecular dynamics in which Ag„ clusters (n = 20- 
200) impacted with total energies in the range of 0.5-6 keV into graphite^i. The penetration 
depth Pd was related to the implantation energy i?o, the cluster number of atoms n, and the 
substrate cohesive energy U according to the following scaling law—: 

oc ^ (1) 

The underlying assumption is that that Eq is mainly spent in breaking bonds of the sub- 
strate while penetrating into it. Therefore, the penetration depth turns out to be inversely 
proportional to the cohesive energy U of the substrate. 

The SCBI technique has been recently applied using Au-nanoclusters (Au-nc) implanted 
on a PDMS substrate^. The resulting nanocomposite can withstand many deformation 
cycles, decreasing their electrical resistance upon cyclical stretching^^. While the effective- 
ness of this novel technique has been demonstratedr^, the physical and chemical phenomena 
underlying the SCBI process at the microscopic scale remains to be fully understood. As 
a matter of fact, the implantation process of neutral metal clusters in a polymer matrix 
is a problem rarely studied in literature^^'i^'SS and the mechanical and electrical properties 
of the resulting nanocomposites need to be fully characterized at the microscopic scale. In 
particular, it is important to investigate the polymer microstructure evolution in terms of 
temperature and pressure waves generated upon implantation and the resulting polymer 
surface damage. Another important property to characterize is the dependence of the clus- 
ter penetration depth on the implantation energy and cluster dimension. This information 



is relevant to better control the SCBI process since it allows to tune the penetration depth 
by changing the implantation energy and/or the cluster dimensions. Finally, it is important 
to assess the different polymer response upon the implantation depending on its microscopic 
structure, namely depending on the entangled-melt vs. cross-linked bond network. 

In this perspective, computer simulations represent a key tool in order to fully describe at 
the atomistic scale both the SCBI process and the following microstructure evolution of the 
substrate. In this work we present molecular dynamics (MD) simulations of the implantation 
of a single Au cluster on a PDMS substrate. We understand that this single event (although 
it can hardly be compared to the full SCBI experiment) is the most fundamental process 
ruling SCBI and, therefore, we believe it deserves a careful investigation. To the best of out 
knowledge a direct simulation of such a cluster impact is still missing in Literature. 

In order to keep the present simulations on the actual experimental length-scale we make 
use of very large-scale MD simulations by aging as many as ~ 4.6 millions atomic trajectories. 
Furthermore, special care is played to carefully model the cluster-substrate interactions, as 
well as the possible presence of molecular linkers into the PDMS film. To this aim, we make 
use of a state of the art force-field which has been accurately benchmarked on well-known 
physical properties of pristine PDMS. 

The goal of the present work is three- fold: (i) modeling at the atomic scale the impact 
(and following penetration) of the cluster into the substrate; (ii) linking the penetration 
depth to the implantation parameters (like e.g. the cluster implantation energy or size, and 
the structure of the substrate); (iii) characterize the surface damage upon impact. 
The much more challenging full simulation of the SCBI process (i.e. the simulation of a 
multi-impact cluster deposition) will be investigated in a following work. 

II. THEORETICAL FRAMEWORK 
A. The force field 

All simulations are performed using the COMPASS (condensed-phase optimized molec- 
ular potentials for atomistic simulation studies) force-field^, including off-diagonal cross- 
coupling terms and high-order (cubic and quartic) force constants. The parametrization of 
the COMPASS force-field allows to describe, using the same functional form, organic (single 
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molecules and polymers) as well as inorganic materials^ii^^. The functional forms used in 
this force field are the same as those used in consistent force-fields (CFF)^: 



total = Eh + Eg + + + Eh^b' + Eh^Q + Eh^^ + Eg^^ + Eg^g/ + Eg^gi^^ + Eq + -EvdW (2) 



where: 
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The functions can be divided into two types: valence terms including diagonal {Ehbi, Eg g') 
and off-diagonal cross-coupling terms {Eh^g, Ei,^^, Eg^^, Eg^g'^^), and nonbond interaction 
terms. The valence terms represent internal coordinates of bond (b), angle (9), torsion angle 
(0), and out-of-plane angle (x), while the cross-couphng terms include combinations of two 



or three internal coordinates (see[T]). The cross-couphng terms are important for predicting 
vibration frequencies and structural variations associated with conformational changes. In 
[3] the subscript denotes the reference values of the bond, angle, dihedral angle and out- 
of-plane angle, while k, ki, k2, fcs, k^ are the force constants determined from the quantum 
mechanical energy surface. Terms involving explicit internuclear distances, r, represent 
nonbond interactions which are composed by a E^dw Lennard Jones (LJ) 9-6 function for 
the van der Waals (vdW) term and a Eg Coulombic function for an electrostatic interaction. 
The LJ-9-6 function is particularly suitable for the case of cluster impacts with respect to 
LJ-12-6 function, which is known to be too hard in the repulsion region. The COMPASS 



force field parameters have been optimized for PDMS by fitting to ab initio calculations2Zi2^. 

B. Entangled-melt PDMS 

We generate amorphous (a) EM-PDMS by means of the Theodorou e Suter— algorithm 
implemented in the Materials Studic^^ package. The simulation cell is composed by 6 chains 
with 100 PDMS monomers each, having in total 6462 atoms and an initial simulation cell 
volume is set as to correspond to a density of 1 g/cm^. After the EM-PDMS structure 
is generated, we equilibrate the cell volume using the LAMMPS package by performing 
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FIG. 1. Graphic illustration of the COMPASS force-filed terms. 
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FIG. 2. Left: EM-PDMS density evolution during the 3.5ns-long NPT equilibration. Right: 
Calculated final radial distribution function. Arrows mark the coordination peaks obtained by 
experimental X-ray diffractio n^^ 



a constant-temperature constant pressure (NPT) simulation for 3.5 ns at a pressure of 1 
Atm and a temperature of 300K. The final density of 0.96 g/cm^ (reported in|2] (left)) is in 
good agreement with the experimental data for PDMS at 300 K— i^i^ii^^. We calculate the 
corresponding radial distribution function g{r) (see|2l right) finding a quite good agreement 
with available experimental datai^"^.— 



C. Crossed- linked PDMS 

The cumputer generation of cross-linked polymers is computationally challenging due 
to the need of simulating the chemical reactions involving the creation of covalent bonds 
between the cross-linker molecule and PDMS. For this reason only few attempts have been 
made in order to describe the PDMS cross-links reactions via MD simulations^"— In this 
work the tetrakis(dimethylsiloxy) silane was used as cross linker (|3]) similarly to Ref.— . 
Tetrakis(dimethylsiloxy) silane is composed by 1 central and 4 peripheral Silicon atoms, 4 
Oxygen atoms, 8 Carbon atoms and 24 Hydrogen atoms. In the presence of a catalyst, 
each tetrakis(dimethylsiloxy)silane cross-linker is able to bond to four terminal hydroxy 
groups on the PDMS chains via its 4 peripheral Si atoms. We simulate this mechanism 
by allowing the covalent bonds between pairs of atoms to dynamically create, during the 
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FIG. 3. Balls and stick representation of the tetrakis(dimethylsiloxy)silane molecule. Carbon 
atoms are shown in cyan, Oxygen atoms in red and Silicon atoms in yellow. The hydrogen atoms 
have been removed for clarity. 



MD simulation, according to specified criteria. In particular, the cross-linking reaction is 
mimicked by forming a bond between one of the 4 terminal peripheral Si on the cross-linker 
and the nearest terminal O on the PDMS chain. The bond forms with unit probability 
when the Si and O atoms are separated by a distance less than 0.1 nm. According to the 
experimental data^-, we use a 8 % cross-linkers weight occurrence, corresponding to 7200 
cross-linker molecules. In order to reduce the computational cost of the PDMS cross linking 
procedure, we initially work on the previously equilibrated EM-PDMS building-block (see 
previous section) by adding 10 cross-linker molecules (corresponding to a 8 % occurrence). 
After this preparatory step, we replicate the so obtained building-block (see below) in order 
to generate the whole PDMS substrate. At this point we performed a NVT simulation 
at T=500 K for 0.1 ns in which the cross-linker molecules and PDMS were allowed to 
dynamically bond according to the criteria specified above. After 0.1 ns more than 90% of 
cross-linkers were bonded to PDMS. 

All the simulations are performed by using the LAMMPS code^. The velocity-Verlet al- 
gorithm with a time step of 1.0 fs is used to solve the equations of motion. A particle-mesh 
Ewald algorithm is used for the long-range electrostatic forces, and the Van der Waals inter- 

8 



actions were cut off at 0.1 nm. The Nose -Hoover thermostat /barostat with corresponding 
relaxation times equal to 100 and 500 fs are used to control the simulation temperature and 
pressure, respectively. 

D. Metal cluster 

As for the description of the Au-nc, we also use the COMPASS force-field, and in particu- 
lar the 9 — 6 Lennard- Jones potential. In fact face-centered cubic (fee) elemental metals can 
be successfully described by means of suitable Lennard- Jones potentials'^. The COMPASS 
9 — 6 Lennard- Jones e and tq parameters for Au have been recently further optimized in 
order to describe: densities, surface tensions, interface properties with organic molecules. 
All these properties are found in good agreement with experiments performed under ambient 
conditions'^. 

E. Simulation procedure 

Once the PDMS building-block (both EM and CL) is equilibrated we obtain the PDMS 
substrate for the SCBl simulation by the following procedure: 

1. We replicate N x Lx M times (where N, L, and M are integers) the basic building-block 
along the x, y and z directions respectively. We then equilibrate such a PDMS bulk 
sample for 0.2 ns at a temperature T=500 K under orthorhombic periodic boundary 
conditions along the x, y and z directions (see H]). This annealing step is crucially 
important in order to remove possible artifacts deriving from the previous fictitious 
periodicity imposed by the rephca of the PDMS basic building-block. 

2. Once the bulk PDMS sample is equilibrated, we create the PDMS surface by removing 
the periodic boundary conditions along the z direction and by further equilibrating 
the slab sample for 0.1 ns at T=300K . While the topmost surface is let to freely age 
under any following MD run, the bottom surface is kept fixed in order to emulate a 
macroscopically thick experimental sample. 



9 



F. Simulation cell architecture 



One of the main issues related to implantation simulations is due to the presence of the 
(x,y) periodic boundary conditions that can eventually give rise to unphysical interactions 
(i.e. overlap between pressure and temperature waves generated upon the cluster impact) 
between neighbor cell replica. In order to avoid such artifacts, we create a frame region 
at the cell boundaries which are coupled to a Nose -Hoover thermostat. The role of the 
frame region is to mimic a bulk absorbing the excess of heat generated upon the cluster 
impact. The rest of the simulation cell, where the cluster impact actually takes place, is 
simulated without any thermostat or barostat. In other words, atomic trajectories in this 
region are aged according to pure Newton dynamics. The simulation cell architecture is 
shown in H] where the fixed bottom region is shown in blue, the frame region coupled to a 
thermostat (lateral and bottom region) in red, while the rest of the cell evolves by means of 
newtonian dynamics. 



TOP VIEW SIDE Vli:W 




FIG. 4. Top (left) and side (right) view of the simulation cell. The fixed bottom region is shown 
in blue, the frame region coupled to a thermostat is shown in red, while the rest of the cell evolves 
by means of newtonian dynamics. The cluster is shown in yellow. 

We then place the Au-nc with the center of mass at a distance of 10 nm from the EM- 
PDMS surface. According to the SCBI experimental data^ the cluster vibrational temper- 
ature before the impact is ~100 K while the PDMS temperature is ~300 K. Therefore we 
equilibrate the Au-nc at T=100K and the PDMS at 300 K for 0.1 ns before the impact. 
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G. Choice of the substrate dimensions. 

We select the optimal substrate dimensions by performing several implantation simula- 
tions on different EM-PDMS substrates of increasing thickness and section. The implanta- 
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FIG. 5. Time evolution of the Au-nc penetration depth for 3 EM-PDMS substrates of increasing 
dimensions corresponding to samples Aq, Bq and Cq. Contrarily to ^-nd Bq, the Cq substrate 
has an optimal lateral dimension and depth that allows to stabilize the cluster position inside the 
EM-PDMS after 180 ps 

tion simulations are performed by considering Au-nc having a radius of 3 nm implanted at 
an energy of 2 eV/atom along the z direction (see H]) . E] shows the time evolution of the 
Au-nc penetration depth during the implantation. The cluster penetration depth at time t 
is calculated as the difference between the maximum EM-PDMS surface height at t=0 ps 
and the cluster center of mass height at time t. We test three different substrate dimensions 
having different lateral dimensions and depths as shown in [H We observed that while the 
absolute penetration depth is not largely affected by the substrate dimension, the long-time 
microstructure evolution sizably depends on them. In particular, while in the case of Aq 
and Bq the cluster position on the PDMS polymer is not stabilized after 180 ps, the Co 
substrate has an optimal lateral dimension and depth that allows to efficiently absorb the 
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Sample Aq 


Sample Bg 


Sample Cq 


x,y (nm) 
z (nm) 
# atoms 


17.07 X nm 

85.35 

1923840 


17.07 X 17.07 
128.01 
2885760 


25.60 X 25.60 
85.35 
4328640 



TABLE I. Lateral dimension, depth and total number of atoms of the substrates Aq, Bq and Cq. 



pressure wave generated upon the cluster impact stabiUzing the cluster position inside the 
EM-PDMS after 180 ps. For this reason we choose the Cq substrate for the implantation 
simulations. The time t = 180 where all the transient post-implantation phenomena are 
terminated will be hereafter identified as tend- 
Ill. RESULTS: IMPLANTATION OF 3 NM AU CLUSTERS 
A. Entangled melt vs. cross-linked PDMS 

After the equilibration of the PDMS+Au-nc system, we perform three different SCBI 
simulation, respectively with implantation energy 0.5, 1.0 and 2.0 eV/atom. We chose Au- 
nc having a radius of 3 nm, comparable to the experimental SCBI data (3 — 10 nm)^^. As far 
as concern the implantation energy, besides 0.5 eV/atom, which is the actual experimental 
value, we choose also larger implantation energies (1.0 and 2.0 eV/atom) to verify the effect 
of this parameter on the cluster penetration and PDMS response. 

[6] shows the time-evolution of the cluster center of mass penetration depth (top) and 
the corresponding velocity (bottom) during the simulation, on the EM-(left) and CL-(right) 
PDMS. 

Based on [6] we can distinguish 4 different microstructure evolution regimes for both EM- 
and CL-PDMS (hereafter referred to as 5*1, 5*2, S'3 and ^4 respectively) depending on the 
polymer response to the cluster implantation. During step 5*1 the cluster moves in vacuo 
by linear uniform motion having a constant velocity depending on the implantation kinetic 
energy. During step 5*2 the cluster impacts on the PDMS surface, penetrates inside the 
substrate and is eventually stopped by friction. During step 5*3 we observe the PDMS 
response to the cluster impact, at first (cluster penetration) the polymer film undergoes a 
sizable compression which is then released by its elastic response, able to drag the implanted 
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FIG. 6. Left: position vs. time of the center of mass (top) and corresponding velocity (bottom) 
on the EM-(left) and CL-(right) PDMS for 3 different implantation energies. 

cluster upward in the polymer. Finally, during step 5*4 the polymer elastic response is 
damped and the cluster stabilizes inside the PDMS substrate at a given depth. These 
four steps are observed for both EM- and CL-PDMS, independently of the implantation 
energy, while the respective time duration every regime takes place depends on the initial 
implantation energy. In particular, we observe that steps 5*2 and are longer by increasing 
the implantation energy, while the cluster stabilization period increases by decreasing the 
implantation energy. 

From|6](top panel) we highlight the differences between the post-implantation evolution in 
EM- and CL-PDMS. We notice that the cluster penetration depth in EM-PDMS is greater 
than in the CL-PDMS. This is due to the fact that in the previous system the polymer 
chains are bonded only via dispersion and electrostatic interactions, while in the latter one 
the chains are partially covalently linked via the cross-linker molecule. Since the cluster 
penetration inside the PDMS matrix involves the breaking of the inter-chains bonds, we 
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conclude that during the penetration inside the EM-PDMS substrate, the cluster experiences 
a weaker friction with respect to CL-PDMS. This picture is confirmed by the analysis of the 
cluster velocity during the implantation (El bottom panel). In the case of the CL-PDMS 
the cluster velocity takes a zero value earlier (~ 35 ps) than in the case of EM-PDMS 
(~ 40 ps). We also observe that the polymer response upon the cluster implantation is 



0.5 eV/atom 1.0 eV/atom 2.0 eV/atom 




FIG. 7. Final structures corresponding to tend of the SCBI simulations with implantation energies 
of 0.5 eV/atom (left), 1.0 eV/atom (center) and 2.0 eV/atom (right). The Au cluster atoms are 
shown in yellow while the PDMS substrate atoms in green. Only the topmost layer of PDMS is 
shown. 

more pronounced in the case of CL-PDMS. In fact, we observe a polymer swelling upon the 
cluster implantation which we quantified as the difference between the polymer maximum 
height after the implantation at t=tend and before the implantation at t = ps. In the 
case of the CL-PDMS we observe a polymer swelling of ~ 3 nm, independently of the 
implantation energy. On the other hand the EM-PDMS swelling was always less than 1 nm. 
This difference can be explained by considering that in the case of CL-PDMS the number 
of internal bonds increases by ~ 8% due to the presence of the cross linker. Therefore, we 
expect a comparatively larger elastic response upon the implantation.lt is worth remarking 
that such a swelling phenomena have been indeed observed experimentally^. This qualitative 
agreement between real and simulated SCBI process stands for the overall reliability of the 
present computational setup. 

[7] shows the final structures at tend of a (x,y) EM-PDMS section having depth of 2 
nm and height of 30 nm (similar structures were observed for CL-PDMS). This cross- view 
confirms that the cluster penetration depth increases with the implantation energy, but also 
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FIG. 8. 3nm cluster penetration depth as a function of the implantation energy (red crosses) fitted 
with a line (green) having angular coefficient of 7 and 6 nm/eV for EM-(left) and CL- (right) 
PDMS respectively. 

adds valuable new information, namely: the polymer surface is strongly damaged upon the 
cluster implantation (in particular for the implantation energy of 2 eV/atom). We observe 
the presence of a crater created on the PDMS surface upon the cluster implantation. The 
lateral dimensions of the crater are almost unchanged by increasing the implantation energy, 
while its depth strongly depends on the implantation energy as well as the surface roughness. 
This feature will be more extensively investigated in a following paper—. 

In [8] we plot the cluster penetration depth vs. the implantation energy in the case of 
EM-(left) and CL-(right) PDMS. As predicted by[T], a linear relationship is found which we 
easily fit finding that by increasing the implantation energy by 1 eV/atom cluster penetration 
depth increases by 7 nm and 6 nm for EM- and CL-PDMS, respectively. This information 
is useful for the SCBI experiments since it allows to tune the cluster penetration depth 
by changing the cluster implantation energy. We further remark that the internal cohesive 
energy contribution Ucl in the case of CL-PDMS is 1.2 times larger than EM-PDMS. There- 
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FIG. 9. Temperature time evolution of the EM- and CL-PDMS for the implantation energies of 
0.5 eV/atom, 1.0 eV/atom and 2.0 eV/atom. The first, second and third lines represents the 
temperature evolutions inside the A, B and C slabs respectively. 

fore we expect a corresponding decrease in in the cluster penetration depth as indeed 
referred in [H 

Besides the cluster penetration depth it is important to characterize the temperature and 
pressure waves generated inside the PDMS substrate upon the cluster implantation. For 
this reason we calculate the substrate temperature during the implantation as a function 
of the PDMS depth. |9] shows the temperature of three slabs {A,B,C) for the EM- AND 
CL-PDMS samples during the 180 ps -long simulation. We choose slab A close to the EM- 
PDMS surface(70-20 nm depth), slab B corresponding the middle of the substrate (370-320 
nm depth), and slab C close to the PDMS bottom (670-620 nm depth). We observe a 
sudden temperature increase upon the implantation for both EM- and CL-PDMS samples. 
In particular the overall temperature increase corresponding to the slab A close to the surface 
(AT4) goes from ~20 K to ~50 K by increasing the implantation energy. This temperature 
wave propagates inside the substrate being strongly dumped for deeper slabs {ATb ~10 K 
and ATq ~5 K). At tend "we observe two main temperature regions present in the polymer, 
namely: a "hot" region (H) at temperature in the range ~ 320-350 K close to the surface 
and another "cold" region at temperature in the range of ~ 300-310 K at larger depth. We 
observe that the overall temperature increase is always larger for the cross-linked substrate. 
We attribute this feature the different stiffness: the deformation upon the impact is larger 
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in the case of the EM-PDMS substrate thus resulting in a comparatively reduced increase 
of temperature. 



t=0 ps t=20 ps 1=36 ps 




FIG. 10. Time evolution of the temperature wave on CL-PDMS during the implantation of a at 
2.0 eV/atom. Each atom is colored according to its local temperature: red color corresponds to a 
temperature of 600 K while blue color corresponds to 300K. The cluster is not shown for sake of 
clarity. 

[TOl shows the time evolution of the temperature wave on CL-PDMS for implantation 
energy of 2.0 eV/atom. For sake of clarity we represent only a PDMS slice having thickness 
of 2 nm and height of 30 nm, without the cluster. Each atom is colored according to its 
local kinetic temperature: red color corresponds to a temperature of 600 K while blue color 
corresponds to 300 K. As due to the cluster impact (for t ~ 20 ps) we observe that the local 
temperature of the atoms in contact with the implanted Au-nc increases up to ~600 K. 
This temperature wave propagates inside the polymer being increasingly dumped with the 
polymer depth. Therefore, besides the temperature increase observed on the PDMS surface, 
we observe a dramatic temperature rising on the cluster impact area. We expect that the 
local properties of the polymer (such as density, number of covalent bonds, pressure) will be 
remarkably changed upon the cluster impact, favoring the penetration of the next clusters 
to be implanted. 

We also calculate the cluster temperature during implantation which suddenly increases 
after the impact as reported in [TTJ In [TT] we summarize the temperature increase ATcz at 
tend which dcpcuds on the implantation energy. The most interesting observation is that the 
structure of the substrate affects the maximum temperature reached at the cluster. 

The cluster implantation gives rise as well to a pressure wave inside the PDMS. [T2] shows 
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FIG. 11. Time evolution of the cluster temperature during SCBI simulation for the EM- (left) and 
CL-(right) PDMS for the implantation energies of 0.5 eV/atom (top), 1.0 eV/atom (middle) and 
2.0 eV/atom (bottom). 





0.5 eV/atom 


1.0 eV/atom 


2.0 eV/atom 


EM - AT^i (K) 


257 


293 


369 


CL - AT,/ (K) 


264 


307 


376 



TABLE II. Cluster temperature increase ATd at the end of the implantation on EM- and CL-PDMS 
samples. 

the component of the stress tensor z [azz along the implantation direction of three slabs 
{A,B,C ) during the 180 ps -long the simulation for the EM-PDMS sample. We choose 
slab A close to the EM-PDMS surface(70-20 nm depth), slab B corresponding the middle 
of the substrate (370-320 nm depth), and slab C close to the PDMS bottom (670-620 nm 
depth). As found for the temperature we observe a pressure wave traveling inside the PDMS 
substrate upon the cluster implantation. From [T2] we can also identify the same — S4, 
steps previously detected: in the first part of the simulation (Si) before the cluster impact 
the pressure oscillates around a constant value. Then, upon the impact (5*2) we observe a 
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Time (ps^ 

FIG. 12. Time evolution of the stress tensor component z (azz) along the implantation direction 
during SCBI simulation for the EM-PDMS. Right panel: the first, second and third lines represents 
the azz evolutions inside the A, B and C slabs respectively. Similar results were found for CL- 
PDMS 

pressure increase on the direction of the cluster initial velocity. Then, we observe the PDMS 
response {S^) marked by a pressure increase on the opposite direction with respect to the 
cluster initial velocity. Finally we observe the regime when the the pressure stabilizes 
around the initial value. Similar results were found for CL-PDMS. 

B. Conclusions 

In the present investigation we study by means of atomistic model potential molecular 
dynamics the implantation of Au-nanoclusters on a PDMS polymer and the following poly- 
mer microstructure evolution. We compare the response of two different PDMS amorphous 
structures: i) entangled melt and ii) crossed-linked. We show that for both entangled-melt 
and cross-linked Polydimethylsiloxane the cluster penetration depth inside the polymer lin- 
early depends on the implantation energy having an angular coefficient of 7 nm/eV in the 
case of entangled-melt PDMS and of 6 nm/eV for cross-linked PDMS. In particular we ob- 
serve that in the case of entangled-melt Polydimethylsiloxane the penetration depth is larger 
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than cross-linked. We show that even a single cluster impact on the Polydimethylsiloxane 
substrate remarkably changes the polymer local temperature and pressure. Moreover we 
observe the presence of craters created on the polymer surface having lateral dimensions 

comparable to the cluster radius and depths strongly dependent on the implantation energy. 
Finally, present simulations suggest that the substrate morphology is largely affected by the 
cluster impact and that most-likely such modifications favor the the penetration of the next 
impinging clusters. 
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